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Abstract 

In the large N c limit, the nonleptonic width of the B meson is determined 
by the lepton mass spectrum in semileptonic B decay and the hadron mass 
spectrum in r decay. This result can be compared with the nonleptonic width 
computed using the operator product expansion and local duality, to estimate 
the violation of local duality in B decay. In the absence of the required data, 
we use experimental Monte Carlo predictions and find that the difference is 
small, at the few percent level. We estimate the theoretical uncertainty in 
our results from higher order corrections (including 1/N C effects). We also 
identify a new contribution to the Aj, — B lifetime difference. 
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I. INTRODUCTION 



Inclusive decays for heavy-quark hadrons such as B mesons are computed in a systematic 
expansion in inverse powers of the heavy quark mass using the operator product expansion 
(OPE). Decay rates are related to the discontinuity of scattering amplitudes at physical 
cuts, so the calculation is necessarily performed for time-like momentum transfer and is not 
a short-distance process for which the OPE is valid. Nevertheless, one can still use the 
OPE for suitably averaged quantities. The simplest example is R(s), the inclusive rate for 
e + e~ — > hadrons, which is related to the discontinuity in the vacuum polarization amplitude 
n(s) along the real axis for s > 0. The operator product expansion can be applied to II(s) 
far away from the physical cuts, such as at complex s = so + £A, A Aqcd, to give H(s) 
as an expansion in Aq CD /s and a s . Analyticity can then be used to relate Il(so + £A) to an 
integral of the form |J 

is. (1) 



:» - so) 2 + A= 



Thus smeared averages of R(s) can be computed using perturbation theory. The a s correc- 
tions to these depend on the smearing width A, and are not under control as A — > 0. In 
practice A is made as small as possible while keeping the radiative corrections under con- 
trol. Typically A is chosen to be A > 500 MeV. One expects that perturbative computations 
agree with data, provided one averages over a large number of hadronic resonances.^ 

A similar situation exists in inclusive semileptonic B decays, B — > X c e v HQ . The dif- 
ferential decay distribution dr < - B,sl - ) /dq 2 dE e dE u can be computed using the operator product 
expansion ||-||. Here q 2 is the invariant mass of the lepton pair, and E e and E v are the 
energy of the electron and neutrino, respectively. The total semileptonic decay rate r( B,sl ) 
automatically involves averaging over different hadronic states, and one expects that the 
OPE computation for r( s ' sl ) is valid. Inclusive nonleptonic B — > X c decays can also be 
computed using an OPE H|J. In this case, there is no variable such as q 2 to integrate over, 
so the computation has no intrinsic smearing. One must then make the assumption of local 
duality, that the quark and hadron computations agree at a single kinematic point. In the 
case of R(s), this would mean that one could compute R(s) for each value of s. This is 
clearly false when there are narrow resonances (e.g., in the T region), but is valid elsewhere 
because the large number of overlapping hadronic resonances provide a natural smearing 
mechanism. An important question is the extent to which local duality holds for nonlep- 
tonic B meson decays. There is a long-standing claimed discrepancy between theory and 
experiment in the B meson semileptonic branching ratio and in the A& to B lifetime ratio, 
and the failure of local duality is one way to resolve these. 

In this paper, we study the violation of local duality using the 1/N C expansion of QCD 
||10|| . In this limit, we show that the nonleptonic B meson decay rate r( B,nl ) can be expressed 
in terms of the product of dr( B ' sl )/dg 2 , the invariant mass distribution of the lepton pair in 



There are examples where local duality is proven to hold even though only two exclusive decay 
modes contribute J2|. 
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semileptonic B decay, and dr( T ' h )/dg 2 , the hadron mass distribution in nonleptonic r decay. 
One gets a quantitative measure of the violation of local duality in the large N c limit by 
comparing the prediction for r( B,nl ) using the OPE with that using the experimental data for 
dT^ B ' s1 ' /dq 2 and dr( T ' h )/dg 2 , and a large N c factorization formula derived here. We estimate 
the accuracy to which the large N c results should be valid in QCD. 



II. BASIC DEFINITIONS AND RESULTS 



Nonleptonic B meson decays are dominantly due to the weak Hamiltonian^ 

c x {m b ) (c llM P L b) (q^P L u) + c 2 (m b ) (q^P L b) (c^P L u)} , (2) 



V2 



aq 



where Pl = (1 — 75)/2 and q = d,s. The Hamiltonian in Eq. (fj) is renormalized at the 
scale /i = m b . The coefficient a( B,nl ) = 1 + (2a em /37r) ln(M w /m b ) is the electromagnetic 
correction. The coefficients q are QCD corrections given by c\ y i = (c + ± c_)/2, where 



C+(fi) 



a s (M } 



-i -6/23 



W j 



C_(/i) 



a s (M } 



-i 12/23 



W i 



(3) 



at one loop order. At next-to-leading order fTTJ, Ci{m b ) = 1.13 and c 2 (m&) = —0.29 for 
a s (Mz) = 0.118. The C\ and c 2 terms are effective charged and neutral current interactions, 
respectively. 

Nonleptonic B decay widths are computed by taking the decay amplitude (X C \H\B), 
squaring, and summing over all possible final states. A simpler method is to compute the 
imaginary part of the B — > B forward scattering amplitude with an insertion of H and W . 
Diagrams involving spectator quarks |fL2,13|, such as Pauli interference, weak annihilation, 
and W-exchange are of order AQ CD /mf relative to free quark decay (since they depend on 
matrix elements of four-quark operators), and can be neglected for the purposes of this 
paper. In the large N c limit, annihilation effects are enhanced by a factor N c relative to 
free quark decay. While this means that annihilation diagrams are formally dominant at 
N c = oo, they can still be neglected for N c = 3 because the N c enhancement does not 
compensate for the 1/mf suppression factor, since the diagram is helicity suppressed and 
involves four-quark operators. 

The free quark decay diagram in Fig. |1] is redrawn in Fig. to show the color flows. 
The charged and neutral current terms in the weak Hamiltonian in Eq. (0) are denoted by 
the exchange of color singlet charged and neutral "gauge bosons." This is a convenient way 
of showing the color flows in the diagrams. The order in N c of the various graphs can be 
computed by using a factor of N c for each closed color loop, and \/y/N~ e for each external 
meson. In the large N c limit, Fig. ||(a) and (b) are 0(N C ), and (c) is 0(1). Fig. ||(a), (b) 



2 The contribution of b — > u decay to the total nonleptonic B decay width is small, and we neglect 
it in this paper. There are also terms in the weak Hamiltonian for b — > cc(d + s) decay, which are 
discussed later. 
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and (c) will be referred to as the charged current, neutral current, and interference diagrams, 
respectively. In the large iV c limit, arbitrary planar gluons must be summed over, which is 
represented schematically by the shaded regions in the figures. The leading diagrams have 
arbitrary planar exchanges within each quark loop, but exchanges between the two quark 
loops are suppressed by at least two powers of N c . Thus in the large N c limit, Fig. 0(a) and 
(b) factorize into the product of two non-perturbative matrix elements. They are: (i) The 
matrix element of two currents in a B meson; (ii) The matrix element of two currents in 
the vacuum. These will be related to inclusive semileptonic B decay and hadronic r decay, 
respectively. 

Consider the inclusive decay of an initial particle / with velocity v due to the Hamiltonian 

The decaying particle I and the weak currents J M and j M will be different for the three cases 
we need in this paper — nonleptonic B decay, semileptonic B decay, and hadronic r decay. 
In these decays, the weak Hamiltonian has factors in addition to those shown explicitly in 
Eq. (|4]), such as CKM angles or renormalization group coefficients, which will be included 
later. 

Define the tensors 

W^(v,q) = (2tt) 3 £<5 4 (m/^ - q - p x ) (I(v)\J^\X) (X\r\I(v)) , (5) 
x 

and 

QT{q) = (2vr) 4 ]T5 4 (g - p Y ) (0\j*\Y) (Y\f\0) . (6) 

Y 

In terms of W and Q, the weak decay rate is given by 



dr 2G 2 F 



dq 2 dq n 



sJq 2 -q*W^(v,q) Q^(q) , (7) 



as long as the diagrams for the process factorize, e.g., as in Fig. 0(a) or (b). This is the 
factorization formula which we will use in the rest of the paper. 

The hadronic tensor W^ v is identical to the one that appears in the OPE for semileptonic 
decays ^ can be expanded as 

w »v = -g^Wi + v*v"W 2 - ie flvaP v a qpW 3 + gVW 4 + {<fv v + qV)W 5 , (8) 

where W{ are functions of q 2 and q-v. Vt^ v can be decomposed into transverse and longitudinal 
parts 

= (-q 2 9^ u + q^q u ) + q^q u , (9) 
where VLt,l are functions of q 2 . Contracting Eqs. (H) and @ gives 

W" % v = Vl T [3q 2 W x - q 2 W 2 + (q ■ vf W 2 

+ Q L \-q 2 Wt + (q- v) 2 W 2 + q 4 W A + 2q 2 (q ■ v) W 5 ] . (10) 
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Both W^ v and Vt^ u are calculable in an operator product expansion. They are related to 
the discontinuities across the cut of the amplitudes T^ u and IL^ defined by 

= -i J d 4 xe~ iqx (B(v)\ T{J^(x) r{0)}\B{v)) , 
W v = -i J d A xe- iC! X (0\T{j^{x) f(0)}\0) , (11) 

i.e., W» u = JmT^/n and Of" = 2lmW / . and IF" have a form factor decomposition 
similar to W^ u and W in Eqs. (|[) and (|]), with W{ — >■ Tj and f2j — > IL> In the complex 
g 2 plane, IF" has a cut along the real axis for q 2 > 0. T^" has two cuts along the real 
axis for q 2 < (m& — m c ) 2 corresponding to intermediate states with a charm quark, and for 
q 2 > (m& + 2m c ) 2 corresponding to intermediate states with a c and two b quarks. 



A. Inclusive hadronic r decay 

The weak Hamiltonian at fi = m T for hadronic r decay r — ► v T X is 

4r y 

# = aW^^^rKsY^w)- (12) 

a( T ) = 1 + (a cm /7r) ln(Mw/ , nv) is the electromagnetic correction. There are no QCD correc- 
tions to Eq. fljjD since the quark part is a conserved current. 

The hadronic r decay amplitude clearly factors into hadronic and leptonic parts, irre- 
spective of whether N c is large or not, since no gluons can be exchanged between the quark 
and lepton parts of the diagram. The leptonic part of the diagram is calculable, and gives 



W[ T) = § (m T - g ) 5(A) , W 2 (r) = m T 5(A) , 
W^ = -W^ = ±5(A), Wf } =0, 



(13) 



where A = m 2 + q 2 — 2m T q . The hadronic part of the r decay diagram is precisely Q for 
the charged current nonleptonic B decay diagram. The factorization formula Eq. (|7|) can 
be applied to this case, with W^ v set equal to Wj$- H in r decay is defined to be fi^L the 
sum of fi's for the weak currents V u d(d , y l j j Pi J u) and V us (s^^Pl u), so that the CKM angle is 
included in tt( T \ Integrating over dgo using 5(A), gives 



dr^ h ) g 2 



2 



(m 2 T -q 2 ) 2 (m 2 + 2g 2 )4 r) +m 2 4 r) , (14) 



dg 2 8ir 2 m%. 

where q 2 is the invariant mass-squared of the hadrons in the final state. 

B. Inclusive semileptonic B decay 

Inclusive semileptonic B — > X c e v decay is due to the weak Hamiltonian 

AC 1 

H = ^f a(B ' Sl) V * (c7^l 6) (eYP L v). (15) 
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a (B,s\) _ j _|_ (a em /7r) ln(Mi^/mfe) is the electromagnetic correction. There are no QCD 
corrections since the quark part of the Hamiltonian is a conserved current. As for hadronic 
r decay, the semileptonic B decay diagram f actor izes into a hadronic and leptonic part, 
irrespective of whether N c is large or not. The leptonic part VL^ U is calculable, 



n 



(B,si) _ 1 (g 2 - m 2 e ) 2 (2g 2 + mj) 



12vr 2q 6 
1 (q 2 - m 2 e ) 2 m 2 e 
8tt q 6 



(16) 



where q 2 is the invariant mass-squared of the lepton pair. Therefore, dr^ s,s ^/dg 2 dgo gives 
information on the hadronic part W^ ,sl \ Using differential B decay spectra for e, /i, and r, 

it is possible to measure w}":*a m. This seems to be a remote possibility, so in the rest of 

this paper we neglect lepton masses. In the m e — > limit, Eqs. (0) and (§) give 



dr^- sl ) _ 2G| 

dq 2 dq 



7T- 



,(B,sl) 



9o - 9 2 



tfW™ + (q 2 - q 2 ) W± 



(B,Bl) 



12tt 



(17) 



2 fi? si) 

where the CKM angle \V c b\ has been absorbed into the definition of Wi 2 ' ■ 



C. The OPE calculation for nonleptonic B decay 



The nonleptonic B decay rate can be computed using perturbation theory if one assumes 
local duality. To order Aq CD /?7?|, the result is 



T^ B ^[b^cu(d + s)] = 3T 



T^[b^cc(d + s)] = 3T c 



,(B.nl) 



2 f C 2 _ + 2c\ 



Ai A | 3A 2 / d__ 2 
2mf J m b ^ dp 



+ 4(c 2 -4)^(l-p) 3 } + 



,(B,nl) 



-> (c 2 _ + 2cl |Y Ax 



2m 2 



_d_ 

d7 



Av 



+ 4(c 2 _ - < ) ^ %) 



mi 



where p = ml/ml, T = G 2 F ml\V cb \ 2 / (192n 3 ) , and 
/(p) = l-8p + 8p 3 -p 4 -12p 2 lnp, 
g(p) = Jl-Ap (1 - Up - 2p 2 - 12p 3 ) + 24p 2 (l - p 2 ) In 



i + y/r^4p- 

1- v/T^4p~ 



/i(p) = Jl - 4p (1 + p/2 + 3p 2 ) - 3p(l - 2p 2 ) In 



i + y/r^4p- 

1- JT^Tp 



(18) 



g(p) 



(19) 



(20) 



The ellipses in Eqs. (0) and JT9D denote corrections of order Ag^/m 3 , and a s {m b ). The 
total nonleptonic decay rate is the sum of Eq. (0) and Eq. (|1^). a s corrections will be 
included when we use these results in the next section. 
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III. A (NOT SO) TOY MODEL: c 2 = 



The total nonleptonic decay rate depends on the short distance QCD corrections C\ and 
C2. The validity of local duality depends on long distance effects, such as the masses and 
widths of hadronic resonances. One can test local duality in the nonleptonic B decay rate 
by studying a hypothetical world in which C2 = and ccid + s) final states are turned off. 
This tests local duality only for the charged current part of Eq. (§). 

The charged current diagram Fig. 0(a) factorizes in the large N c limit, so Eq. (|?p can be 
used. Combining this with Eqs. (|10D, (|I4]), and flTT]), one finds 



dr( B < nl ) 

dq 2 



a 



(B,nl) 



a 



(B,Bl) a (T) 



dq 2 a 



(B,Bl) 

T 



4 r) (<? 2 



dr( T - h ) 96yr 3 m! 



(<? 2 ) 
1 



n 



(B,sl) 



(q 2 



dq 2 



r<2 



{mi 



q 2 ) 2 (m 2 + 2q 2 



(21) 



Ql has been neglected in this equation because it is proportional to the current quark masses 
of the u, d and s quarks. In the large N c limit, one can compute the nonleptonic B decay 
rate in terms of the semileptonic B decay rate and the r hadronic decay rate. This is the 
main result in this paper. 

Equation ([H]) will be used to test local duality in nonleptonic B decay. Assuming 
that the OPE computation of the total semileptonic decay rate is valid, one knows that 
the integral of the experimentally measured spectrum dT^ B ' sl ^/dq 2 over q 2 is equal to the 
OPE value of the total decay rate. Moments of the measured q 2 spectrum and of the OPE 
prediction should also agree; however, the higher order a s and 1/m corrections get larger 
for higher moments. As a result, the integral of the q 2 spectra from experiment and from 
the OPE should agree when weighted with a sufficiently smooth and broad function, but 
not necessarily point by point. A similar result holds for the q 2 spectrum in r hadronic 
decay. In Eq. ( |2l| ) above, the experimentally measured q 2 spectrum in semileptonic B decay 
is multiplied by the experimentally measured q 2 spectrum in hadronic r decay. The r decay 
spectrum is not a smooth function (see Fig. 3), which is what makes the test non-trivial. 



A. Error estimates 

The test of local duality is to compare T^' 111 ) in Eq. ([ID calculated using the OPE 
and perturbation theory with r( B,nl ) in Eq. ([H|) obtained from experimental data. Before 
we discuss the result, it is worth examining the errors in the formulae we have derived. 
The leading corrections to Eq. fl2"T| ) are factorization violating diagrams involving two-gluon 
exchange between the two quark loops in Fig. |2|(a). Single gluon exchange is forbidden by 
color conservation. If both gluons are hard, the diagram can be computed in perturbation 
theory. It is order a 2 , but is not part of the BLM series [16[ {a s /^) r+1 where /3q = 



11 — 2nj/3 is the coefficient of the one-loop QCD /5-function. If both gluons are soft, the 
correction is of order Aq CD /m^, and if one gluon is hard and one is soft, the correction is 
of order a s AQ CD /m^. The leading power suppressed corrections to Eq. (|2lD are of order 
a s AQ CD /m^ or AQ CD /mf. The AQ CD /m^ corrections cancel between r( B,nl ) and r( B ' sl ). The 
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first correction in the OPE for is order AQ CD /g 4 . Both r decay and B decay involve 
summing over du and su final states with weights \V u d\ 2 and IKJ 2 , respectively, so both 
contributions are correctly included in Eq. (|21| ) when one uses the measured r hadronic 
decay rate. In a world with C2 = 0, the largest correction to Eq. (|2l|) is the Pauli interference 
contribution,^ of order (167T 2 / 7V C ) f B / 'm 2 B ~ 5%. Very little use has been made of the 1/N C 
expansion in the theoretical analysis for c<i = 0. However, the skeptic might argue that the 
error estimate is based on the OPE, which is what is being tested in the first place. The 
factorization violation contributions to the decay rate are 1/A^ 2 corrections. This gives an 
OPE-independent error estimate of ~ 10%. 



B. Data analysis 

The hadronic r decay spectrum (i.e., il^) is only measured in the region q 2 < 



m 



2 



In the region A 2 < q 2 < m 2 B , (with A 2 < m 2 ), we use the perturbative calculation ||14 



including corrections to ^ T \q 2 ) up to order a 2 with a s (m T ) = 0.32. Higher order terms are 
negligible for our purposes. We vary A 2 between 2 — 3 GeV 2 and check that our conclusions 
are unaffected by the precise value of A. The region q 2 < m 2 is where the largest violations 
of local duality are expected in Q, and experimental data suggests that local duality already 
holds at some level for q 2 > 1.5 GeV 2 . For the largest deviations from local duality are 
expected near g 2 ^ ~ {m h — m c ) 2 . 

For the hadron invariant mass spectrum in r decay we use a Monte Carlo distribution 
that has been tuned to be consistent with CLEO data.0 This can be used to determine 
Vtj,\q 2 ). In Fig. |3]we have plotted VLj.\q 2 ) /Vl!^ ,sl \q 2 ) together with the perturbation theory 
prediction for this ratio (including corrections up to order a 2 ). The lepton invariant mass 
spectrum in semileptonic B decay has not yet been measured, so we use the CLEO Monte 
Carlo.Q The normalized mass spectrum is plotted in Fig. |as a function of q 2 together with 
the prediction of perturbation theory. The shape of this curve is rather insensitive to the 
order a s and a 2 (3o corrections [fl8| . 



It is simplest to normalize the nonleptonic rate to the semileptonic rate. This eliminates 
the order AQ CD /m 2 corrections and reduces the perturbative corrections as well. The result 
for the B meson decay rate using the factorization formula Eq. (|2lD and the CLEO Monte 
Carlo spectra is 



T^[b^cu{d + s)}/T^ B 



,sl) 



CLEO 

The result for the OPE calculation at this order is 



3.32 |ci| 2 . (22) 



3 It is possible that the non-BLM a 2 correction, which has not yet been computed, is large. 

4 We thank A. Weinstein for the Monte Carlo spectrum. Our results would not change if we used 
the ALEPH measurement [17|. 



'We thank E. Potter and A. Weinstein for the Monte Carlo spectra. 



dr( s ' nl )[&^ cu(d + s)]/dq 2 




2 " 


dr( B ' sl )/dg 2 


a (B,s\) 





7T 



0.173/3 




(23) 



Expressing a s (g 2 ) in terms of « s (m;,) = 0.22 and neglecting nonperturbative and 



higher order perturbative corrections, this implies T^ B ' nl ^[b — ► cu[d + s)]/T^ B 



si) 



C\CL 



(B.nl) / Q (B,sl) 



1 + 0.07 + 0.025) = 3.21 Idf. The difference is surprisingly small, at 
the 3% level, which is less than the 5% theoretical uncertainty of the calculation. So we 
conclude that the uncertainties in the charged current contribution to b hadron lifetimes 
due to local duality violation are unlikely to exceed the ~ 5% percent level. Of course, it 
is very important to repeat this analysis using real data for the lepton mass spectrum in 
semileptonic B decay, dr( B,sl )/dg 2 , especially since the Monte Carlo predictions used above 
(solid curve in Fig. M differ significantly from perturbation theory in the low q 2 region. 



IV. REALITY 



In the real world with c 2 7^ 0, there are several effects that restrict our calculational 
ability. The interference term Fig. 0(c), is a 1/N C correction and does not factorize. Its 
contribution is order (2/iV c )|£^/ci| ~ 18% of the B — > cu(d + s) decay rate. It has been 



argued that deviations from local duality are likely to be larger for this contribution [IS . 

The neutral current decay Fig. |](b) also factorizes. W^ v for this decay can be determined 
from semileptonic B decay due to the (hypothetical) neutral current dj^P^b. This can be 
related to semileptonic B decay due to the weak current vry^Pifi by isospin invariance. 
for the neutral current decay is the vacuum polarization contribution due to the cy^Pjju, 
current, for which there is no experimental information. The neutral current contribution 
to the decay width is order |c 2 /ci| 2 ~ 7%. 

We have neglected b — » cc(s + d) decays. These can also be factorized in the large N c 
limit. The vacuum polarization Q is due to s^P^c and d^P^c currents, for which there is 
no experimental information. Thus our calculational method, while valid in principle, is not 
useful for b — > cc(s + d). Experimentally, one can separate the b — > cu(s + d) and b — > cc(s+d) 
contributions to the nonleptonic decay rate by counting the number of charm quarks in the 
final state, and the results of the previous section can be used for the b — > cu(s + d) part. 



V. CONCLUSIONS 

We have estimated the violation of local duality for nonleptonic B decay with c 2 = 
to be less than the theoretical uncertainty in our calculation of ~ 5%. It is important to 
redo the analysis using experimental data, rather than Monte Carlo, once it is available. 
One can also view the results in a different way: the factorization formula allows one to 
compute the nonleptonic B decay rate using the lepton mass spectrum in semileptonic B 
decay and the hadron mass spectrum in r decay without using an OPE. In the world c 2 = 0, 
one can predict this rate with an accuracy of 5%. In the real world, one can compute the 
b — ► cu(s + d) part of the decay rate with approximately 20% accuracy. The results cannot 
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easily be extended to D meson decays, because the A.Q CD /ml corrections in the charm sector 
are ~ 100% of the charm quark decay result. 

These reasonably tight limits on local duality violation do not resolve the question of 
the Ab to B lifetime ratio. Of course, duality violation in b — > cc(d + s) and in the C1C2 
interference contributions are probably larger. While enhancing the average number of 
charm in B decay seems experimentally disfavored, the c\Ci interference remains suspect. 
It is also possible that spectator effects are anomalously large in decay — in fact, there 
are six-quark operators (see Fig. 5) which contribute roughly of order {IQir^^fs/mbY/Nc 
relative to the b quark decay. These have not been analyzed, and may significantly affect the 
spectator contribution to A b decay. They do not occur for B meson decay, and so contribute 
to the Ab — B lifetime difference. 

In semileptonic B decay, there is no 1/rrit, correction p. This is also true in nonleptonic 
B decays, if one uses the OPE 0. In other processes which do not have an OPE, it is known 
that one can get arbitrary power suppressed corrections [^TJ. It has been suggested that 
there might be 1/rrif, corrections in inclusive nonleptonic decay rates [[H]]. In the large N c 



limit with C2 = 0, Eq. ([H]) shows that the 1/rrib corrections to the nonleptonic decay rate 
can be computed in terms of the corrections to the lepton mass spectrum in semileptonic B 
decay and the hadron mass spectrum in r decay. At least one of these would have to have 
a l/rrify correction if there is one in the nonleptonic decay rate. 
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FIG. 1. Leading contribution to B meson decay (B — > B forward scattering amplitude). The 
weak currents have been contracted to a point, and (g> represent quark bilinears that create or 
annihilate B mesons. 




(c) 

FIG. 2. Free quark decay diagrams showing the color flow. The charged and neutral current 
interactions in Eq. (||) are represented by the exchange of gauge bosons denoted by wiggly and 
dashed lines, respectively. The shaded regions represent arbitrary planar gluons. 
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FIG. 3. Q)p' (q 2 )/Q, T ' (g 2 ) from Monte Carlo fitted to r decay data (solid curve), and the 
prediction of perturbation theory (dashed curve). The large resonance peak is due to the p. 



0.2 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 r 




q 2 (GeV 2 ) 

FIG. 4. Monte Carlo (solid curve) and perturbation theory (dashed curve) predictions for 
(l/r( B ' s1 )) dr( B ' sl )/dg 2 . The area under each curve has been normalized to unity. 

d > u 

b >■ ^ b 

u * d 

FIG. 5. Spectator diagram which contributes to the lifetime, but has no analog for B meson 
decay. 
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